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Shuttling of electrons in single-molecule transistors with magnetic leads in the presence of an 
external magnetic field is considered theoretically. For a current of partially spin-polarized electrons 
a shuttle instability is predicted to occur for a finite interval of external magnetic field strengths. 

The lower critical magnetic field is determined by the degree of spin polarization and it vanishes 
as the spin polarization approaches 100%. The feasibility of detecting magnetic shuttling in a 
Ceo-based molecular transistor with magnetic (Ni) electrodes is discussed [A. N. Pasupathy et al., 

Science 306, 86 (2004)]. 
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1. INTRODUCTION 

In recent years the effect of the spin of electrons on the 
transport properties of nanostructures have been studied 
intensively, both theoretically and experimentally. In the 
context of spin-based electronics (spintronics) the possi¬ 
bility to control electrical currents by a weak external 
magnetic field using the Zeeman and/or the spin-orbit 
interaction is one of the main goals. 

Magnetic materials and especially half-metals are nat¬ 
ural sources of spin-polarized electrons for spintronics. 
Transport of spin-polarized electrons in nanostructures 
(quantum dots, suspended nanowires, etc.) in exter¬ 
nal magnetic field results in new phenomena where spin, 
charge and mechanical degrees of freedom are strongly 
inter-related. In this new field of investigations (spin- 
tromechanics, see Ref. [2) the presence of a mechanically 
“soft” subsystem results both in a strong enhancement of 
spintronic effects and in magnetic control of the mechan¬ 
ical subsystem in the classical as well as in the quantum 
transport regimes. 

Vibrational effects are known to be important for the 
transport properties of molecular transistors (see, e.g., 
the reviews in Refs. 0 and 3). In single-molecule transis¬ 
tors a strong electron-vibron coupling was observed in a 
Cgo-based transistor with nonmagnetic (gold) leads Q- 
The measured current-voltage characteristics in this ex¬ 
periment revealed low-energy periodic step-like features. 
They were interpreted as a signature of vibron-assisted 
electron tunneling via the fullerene molecule. Experi¬ 
mental I — V curves were theoretically explained in 
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the frames of a simple model of a single-level quantum 
dot strongly coupled to a single vibrational mode and 
weakly coupled to the source and drain electrodes. 

Later on Cgo-based molecular transistors with mag¬ 
netic (Ni) leads were fabricated fj;]. In samples where 
the tunneling coupling to the ferromagnetic electrodes 
were relatively strong (~ tens of meV), Kondo-assisted 
tunneling via the Cgo molecule was observed. These mea¬ 
surements also proved the presence of a strong inhomoge¬ 
neous magnetic field produced by the ferromagnetic elec¬ 
trodes in the nano-gap between them. In samples with 
weak tunneling couplings the usual Coulomb blockade 
picture for a single-electron transistor was observed. 

In the present paper we formulate the conditions for 
the appearance of a vibrational instability of a fullerene 
molecule suspended in the gap between two magnetic 
leads with opposite magnetization. Electron shuttling 
of spin-polarized electrons produced by magnetic (ex¬ 
change) forces was predicted in Ref. 0 for the case of 100% 
polarization of the leads. In this limit (realized for half¬ 
metals) the electric current is blocked (spin blockade) 
in the absence of spin-flips induced by, e.g., an external 
magnetic field. It was shown that in the absence of dissi¬ 
pation in the mechanical subsystem such a magnetic field 
triggers a shuttle instability even for vanishingly small 
fields 0. In the presence of dissipation a threshold mag¬ 
netic field is determined by the rate of dissipation and it 
is small for a weak dissipation. 

One of our aims here is to develop a theory of magnetic 
shuttling for conditions corresponding to the experimen¬ 
tal set-up of Ref. where the electrons in the ferromag¬ 
netic leads were partially polarized (~ 30%). The ab¬ 
sence of a spin blockade in this case qualitatively changes 
the criterion for electron shuttling. We will show that 
even in the absence of mechanical dissipation, a shuttling 
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regime of electron transport occurs in a finite interval of 
external magnetic field strengths, H m j n < H < H max , 
where H m ; n is determined by the degree of spin polariza¬ 
tion 77 . 

In particular, for a high degree of spin polarization 
(77 —>■ 1) and if T Huj (where T/H is the tunneling rate 
of majority spin electrons and ui/2ir is the mechanical 
vibration frequency of the fullerene) the threshold mag¬ 
netic field for reaching the shuttling regime of electron 
transport reads: H m \ n ~ ^/l — 77 T. In the limit of “hard” 
vibrons, T -C tux, the threshold field is determined by the 
vibron energy, H min ~ yj 1 — 77 tux. 

The calculations outlined below aim at determining the 
rate of change r(H) (to be defined later) of the center-of- 
mass coordinate of the single-molecule shuttle, the sign of 
which then allows us to formulate the conditions required 
to observe shuttling of spin-polarized electrons in Ceo- 
based molecular transistors. 


2. HAMILTONIAN AND EQUATIONS OF 
MOTION 


The Hamiltonian of a magnetically driven single elec¬ 
tron shuttle (see Refs.0 Gd]) consists of four terms, 

H= +H d + H v , ( 1 ) 

j=S,D 


where Hj is the standard Hamiltonian of noninteract¬ 
ing electrons in the source (j = S ) and drain (j = D) 
electrodes, Htj is a tunneling Hamiltonian with coor¬ 
dinate dependent tunneling amplitudes tj (x) and H d is 
the Hamiltonian of a single level (eo) quantum dot (QD) 
magnetically coupled to leads of spin-polarized electrons 
by coordinate dependent exchange interactions Jj(x). It 
has been shown 0 ] that the shuttling regime of single¬ 
electron transport can be realized in the presence of an 
external magnetic field H ext for oppositely magnetized 
source and drain electrodes. If the external magnetic 
field is directed perpendicular to the antiparallel polar¬ 
ization vectors of the leads, the QD Hamiltonian reads 


H d = 


£0 - 

g[iH 


J(x) 




£0 


J(x) 


a \ a i 


(ajaj. + aja t ) + U aja-|-aja,|,. (2) 


Here a£ (a CT ) is the creation (annihilation) operator for 
an electron with spin projection a =j\ -l on the dot, 
J(x) = J s {x) — Jd(x), p is the Bohr magneton (g is 
the gyromagnetic ratio), and U is the Coulomb repulsion 
energy. QD vibrations are described by the harmonic 
oscillator Hamiltonian 


-9 2 ~ 2 

- P mu x 
H ’=2^ + —- 


(3) 


where x is the displacement operator, p is the canonical 
conjugated momentum {[x,p] = ih), m is the mass and 
u is the (angular) vibration frequency of the QD. 

The aim of the present paper is to find the conditions 
under which magnetically driven shuttling in experiments 
with fullerene-based single-molecule transistors 00 can 
be realized. In this case the ferromagnetic leads are char¬ 
acterized by a certain degree 0 < 77 < 1 of spin polariza¬ 
tion (~ 30% in the experiment of Ref. 0) and the charac¬ 
teristic vibron energy tux (~ several meV, Ref. 0) is larger 
or of the order of the energy scale T that characterizes 
the tunneling coupling to the leads, T < tux. 

It follows from the above considerations that we are 
not in the adiabatic regime of mechanical motion (which 
is said to be antiadiabatic when T <C tux) and the phys¬ 
ical picture of “magnetic shuttling” developed in Ref. 0 
(see also Ref. I 111) fo r the adiabatic regime does not hold 
here. As in Ref. Il2| we will solve the problem using equa¬ 
tions of motion for the reduced density operator. In the 
limit eV T, tux, pH, ksT (T is the temperature, V is 
the bias voltage) the density operator can be factorized 
into the product of a QD density operator and an equi¬ 
librium density matrix for the leads. In the Coulomb 
blockade regime, eV. T <C U, the matrix elements of 
the QD density operator (p 0 = (0|/5 d |0), p a = (ojp d |cx), 
PIT = (t \p d \ 4-)) in the Hilbert space of a singly occupied 
dot level (they are still operators in the Hilbert space of a 
harmonic oscillator) are determined by the following set 
of equations 0, 

= -i[ H v,Po]~ \ {r1;(x) + T^(x),po| 

+ \/ r L(£)PT\/ r L( i ) + \J T D^)Pi\J T D^)^ ( 4 ) 

= -* [ H v,Pt\ + \ Pt\ + y (m - Pn) 

+ y/Tl(x)po^Tl(£) - i {r£(x),pt} > (5) 

-§j~ = ~ i l H v’Pi\ - \ [ J ( x )>Pi\ - y (rn - p\i) 
+\l T s( x )Po^/ T s( x ) ~ \ {• . ( 6 ) 

yy 1 = -i [H v ,p n ] + l - {J{x),p n } + y (pt - Pi) 
-yM,(z)- (7) 

In order to solve Eqs. Q 0 it is convenient to in¬ 
troduce dimensionless variables for time (tix —> t ), dot 
displacement ( x/xq -a x, where Xq = yjti/mix is the 
zero-point oscillation amplitude), momentum ( pxo/h —> 
p ) and for various characteristic energies (tux —> 1, 
g pH/tux —> h, TJ(x)/tux -A TJ(x), where TJ(i) = 
2m/\tj >a (x)\ 2 is the level width, v is the density of states). 
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3. INFLUENCE OF SPIN POLARIZATION ON 
THE MAGNETIC SHUTTLE INSTABILITY 

We are interested in the conditions under which a mag¬ 
netic shuttle instability occurs. Therefore it is sufficient 
to expand tunneling amplitudes and exchange energies 
up to linear terms in the coordinate operator x, 

J(x) = Jq — ax, Tg/ D {x) = r a S / D T j (8) 


for which ImO < 0 , corresponding to the shuttle regime 
of electron transport, satisfy the inequality 

h 4 + A x h 2 + A 2 < 0, (13) 


where 

Ai 

A 2 


-3 (rf +1) - R_It (sr? + 4), 

(r?-rl)(rf + 4) / r?-r|\ 

8Tf + 1 + 


(14) 

.(15) 


where l is the characteristic electron tunneling length (of 
the order of 1 A). 

Note that the zero-point fluctuation amplitude of a 
fullerene molecule in the vicinity of the Lennard-Jones 
potential minimum is much smaller than the electron 
tunneling length, l 1. Since the spatial scale of the 
exchange interaction is also determined by the tunnel¬ 
ing length l, we can treat the mechanical coordinate x 
as a classical variable. The equation of motion for the 
classical coordinate x c takes the form 


d 2 x c 

dt 2 


+ x c 


-| Tr (Pt-Pl) • 


(9) 


The r.h.s. of Eq. © can readily be found from the set of 
linear equations 0 0 when the operator x is replaced 

by x c . 

To simplify analytical calculations we consider a sym¬ 
metric junction for which Jo = 0 , Tg = T^ = T, and 
Ts = rjj = 7 . To first order in the small dimensionless 
parameter a/l ~ ( J/Huj) (xq/1) 2 <C 1 Eq. 0 is reduced 
to a linear equation. For the time dependent coordinate 
displacement x\ = x\ (t) this equation reads 


d 2 x i 
dt 2 


+ Il = ~2l 



dt'(eo|e ^ 4 t ^\e)xi{t'). 


( 10 ) 


Note that in the absence of an external magnetic field 
the nanoelectromechanical coupling results in additional 
damping of the mechanical subsystem (since the coeffi¬ 
cient A 2 > 0 ). 

It is evident from Eq. m, which is biquadratic in the 
magnetic field ( h ), that a shuttle instability occurs in a 
finite interval of magnetic fields, h m i n < h < /i max . Now 
we introduce the degree of spin polarization as 


r-7 

N t + N i r + 7 


(16) 


where N a is the number of particles with spin projection 
<j. The last relation in Eq. (fTfil) comes from the definition 
of tunneling rates and our assumption that tunneling am¬ 
plitudes do not depend on spin projections. We analyze 
Eq. (THU) in two limiting cases: (i) T> 1, 1 — 1 , and 
(ii) T « 1 , 1 — 77 « 1 . 

Note that in the limit of weak polarization, rj —) 0, 
magnetic forces are small, J —> 0, and electron shuttling 
is supported by Coulomb forces. This case was consid¬ 
ered in Ref. [l2|, where it was shown that in the limit 
of weak polarization a magnetic field ceases to influence 
single electron shuttling. 

In the adiabatic limit, (i), one finds (reverting to using 
parameters with dimensions) that the shuttle instability 
region is defined by the double inequality 


Here the matrix A takes the form 



( 11 ) 


Ti^ = T ± 7 , |eo) = ( 1 , 0 , 0 ) T , and the vector |e) = 
(ei,e 2 ,e 3 ) T has the components 


<Vsr. (17) 

In the antiadiabatic limit, (ii), the critical magnetic fields 
are determined by the vibron energy hio rather than the 
tunneling rate T, so that 


hio 



< ggH < VZ fuo. 


(18) 


= -fu; - Ii), = -fi (r? - + 4)7 

e 3 =- ^-=-, A = 3T 2 + r 2 + 12h 2 . (12) 

A shuttle instability occurs when the oscillatory solu¬ 
tion of Eq. m, x\ (t) ~ exp(iflt) becomes unstable, that 
is when Imfl < 0. The smallness of the r.h.s. of Eq. m 
allows one to obtain the imaginary part of fl by pertur¬ 
bation theory. The values of the parameters Ti, T 2 and h 


For the general case that part of the T,H parameter 
space which corresponds to a magnetic shuttle instabil¬ 
ity is shown in Fig. 1 for different degrees of spin po¬ 
larization. For a given (normalized) tunneling coupling 
T /huj the range of magnetic fields for which an instability 
occurs is shifted towards higher fields as the spin polar¬ 
ization decreases. The upper critical field for high spin 
polarizations, rj —> 1 , depends linearly on rj , 

ff^max ~ y/3 (T 2 + (M 2 ) [1 + K(1 - g)} , (19) 
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FIG. 1: Upper critical magnetic field (upper set of curves) 
and lower critical field (lower set of curves) for the shuttle 
transport regime plotted as a function of the normalized tun¬ 
neling rate T/hw of majority spin electrons for different values 
of spin polarization g. The solid curves were plotted for g = 1 
[100% spin polarization, here the lower curve coincides with 
the a;-axis], the dashed curves were plotted for g = 0.8, and 
the short-dashed curves for g = 0.3. The shuttle regime corre¬ 
sponds to the area between the lower and upper critical fields 
(dark region for the case of g = 0.8), while outside this area 
one is in the vibronic regime. 


where K ~ 1 is a positive constant. 

The dependence of the lower critical magnetic field on 
g is weaker. In the antiadiabatic regime one finds that 

q^H min ~ hwyj 1 - 77 , T < hw, ( 20 ) 

and hence H m j n rapidly saturates to a constant value of 
order hw/gg with decreasing spin polarization (see Fig. 2, 
short-dashed curve). In the adiabatic regime, T hw , 
the lower critical field decreases linearly with increasing 
spin polarization, except in the close vicinity of complete 
spin polarization, where H m j n oc y/1 — g (compare the 
solid curve in Fig. 2). 

The appearance of an upper and a lower critical mag¬ 
netic field has a simple physical explanation. When gH 
is the largest energy scale in our problem, gH F, hw, 
the fast precession of the electron spin of the dot in a 
perpendicular external magnetic field nullifies the aver¬ 
age spin and the magnetic shuttle instability disappears. 
To estimate the upper field one may compare the char¬ 
acteristic spin precession frequency, gH max /h, with the 
electron tunneling rate, T/h, or the frequency of vibra¬ 
tions w. That is gH max ~ max (T, hw). The lower critical 
field can be readily estimated for a high degree of spin po¬ 
larization, 1 — r] -C 1. In this case we have to compare the 
average time between spin flips, r/, induced by a constant 
magnetic field H in the presence of an electron tunneling 
coupling r with the characteristic life-time of minority 



FIG. 2: Plots of the lower critical magnetic field H m i n , which 
defines the border between the vibronic regime (H < ffmin) 
and the shuttle transport regime (H > if min; compare Fig. 1) 
vs. spin polarization g for different values of the normalized 
tunneling rate T /hw of majority spin electrons (solid curve: 
T/hw = 10; dashed curve: T/hw = 3; short-dashed curve: 
T/hw = 0.1). 


spin electrons on the dot, ~ h/'y. The spin-flip rate Vf in 
weak magnetic fields H can be estimated by perturbation 
theory with the result that hvf ~ (nH) 2 /max(T, Hui). 
Therefore the lower magnetic field is strongly sensitive 
to spin polarization, 

~ -\/f 7 max(r, hui) ~ y/l — r] max(r, Hui), (21) 

and disappears for 100 % spin-polarized electrons (77 = 1 ). 

Next we estimate the maximum rate of (exponential) 
increase, r m = — Im{fI(Ff op t)}, of the QD oscillation am¬ 
plitude in the shuttle regime. In the adiabatic limit, 
T » hw, one finds that ggH opt ~ 0.4 F and that 

r„ = C f(5) 2 . (22, 

where C ~ 0.1 is a small numerical factor. In the case 
r <C hw, which we are interested in here, the maximum 
rate is realized when ggH 0p t ~ hw, corresponding to 


~ f ^£.) 2 

h hw V 1 ) 


(23) 


where we omit a numerical factor of the order of one. 

In the presence of dissipation in the mechanical sub¬ 
system, which can be described by adding a phenomeno¬ 
logical friction term 7 d,x c (t) to the equation of motion 
© ('Yd = u/Q, where Q is the quality factor), the shut¬ 
tling regime appears when r m > w/Q. Therefore electron 
shuttling in a Ceo-based molecular transistor with mag¬ 
netic electrodes could be realized if the quality factor Q 
of the mechanical resonator obeys the inequality 


For the experimental setup in Ref. 0, where fullerene 
vibrations were observed, the factor (l/x o ) 2 — 10 3 and 
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r -C Hco ~ 5 meV (one can estimate T ~ 0.1 0.5 meV 

from the maximal current measured in Ref. 5 ). In 
the C 60 -based transistor with magnetic (Ni) leads J ~ 
r ~ 10 meV (see Ref. 0 ). From Eq. EH) one can esti¬ 
mate that the required quality factor is Q > 10 3 — 10 4 . 
However the optimal external magnetic field in this case, 
H opt — 50 T, is too high. Instead, we therefore estimate 
Q for magnetic fields in the vicinity of the lower critical 
magnetic field H > H mm where magnetic fields for a very 
high degree of electron spin polarization (~ 99%) could 
be of the order of a few tesla. In this case (Hui F, 
1 - r] -C 1) 


jr(i - rf) (j_\ 

- r 2 + 4( 1 -7 ? )(M 2 Uoj 


(25) 


Assuming that T ~ tko we find Q ~ Qopt/(l — rf). 


4. CONCLUSIONS 

In summary we have considered the feasibility of ob¬ 
serving magnetically driven single-electron shuttling un¬ 
der realistic conditions corresponding to an already ex¬ 
perimentally realized Cgo-based single-molecule transis¬ 
tor with magnetic leads. The main requirement for mag¬ 
netic shuttling is the presence of an external magnetic 
field that induces electron spin flips. We have shown 
that the optimal magnetic field, defined as the field that 
maximizes the rate of increase of the shuttling ampli¬ 
tude, is determined by the vibration frequency w. For 
fullerene-based single-electron transistors this frequency 
could be in the THz region [ 4 ] with corresponding mag¬ 
netic fields in the region of several tenths of teslas. For 
magnetic electrodes with a very high degree of spin polar¬ 
ization one needs less strong (by an order of magnitude) 
magnetic fields. However, the quality factor of the cor¬ 
responding mechanical resonator has to be exceptionally 
high, Q > 10 5 . 
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